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ABSTRACT 

Using cosmological simulations, we make predictions for the distribution of clusters 
in a plausible non-gaussian model where primordial voids nucleated during inflation act 
together with scale- invariant adiabatic gaussian fluctuations as seeds for the formation 
of large-scale structure. The parameters of the void network are constrained by the 
cosmic microwave background (CMB) fluctuations and by the abundance and size of 
the large empty regions seen in local galaxy redshift surveys. The model may account 
for the excess of CMB temperature anisotropy power measured on cluster scales by 
the Cosmic Background Imager (CBI). We show that the z — cluster mass function 
differs little from predictions for a standard ACDM cosmology with the same <j$, 
but that the evolution of the mass function at z ~ 1 is slower than in a gaussian 
model. Because massive clusters form much earlier in the "void" scenario, we show 
that future integrated number counts of SZ sources and simple statistics of strong 
lensing will provide additional constraints on this non-gaussian model. 
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1 INTRODUCTION 

The emergence of the ACDM scenario as the standard 
model for both the evolution of the cosmological back- 
ground and the development of large-scale structure has re- 
cently received dramatic confirma tion from the WMAP mi- 
crowave background experiment iSpergel et aljj2 003f) and 
from the large-scale distrib ution of galaxieT TPearocket~ai] 
l200ll:lTegmark et aljl2003ef) . 

In this picture of minimal complexity, adiabatic gaus- 
sian fluctuations in the energy density with scale-invariant 
power law power spectrum generated in the early Uni- 
verse are stretched during inflation to astrophysically rel- 
evant scales. Gravitational instability is then responsible 
for amplifying the resulting overdensities and for the col- 
lapse of structure. In hierarchical models like cold dark mat- 
ter (CDM), small objects collapse first and the more mas- 
sive clusters form relatively recently. Although successful in 
many aspects, the ability of the model to re produce the 
mass function of satelli tes of the Milky- Way (iTullv et alJ 
l2002l:IStoehr et alJl20o3h the rotation curves of LSB galax- 
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ies iMcGaueh et al.l2003lh or the bulk properties of t he stel- 
lar populations of massive ellipticals <Peebledl200lh is still 
uncertain, even if recent hints to wards tilt and /or r unning 
of the primordial spectrum index dPeiris et alJl2003h might 
help resolve some of the discrepancy on small scales. 

On larger scale s, the voids seen i n the nearby 
galaxy distribution JEI-Ad fc Piranl 1200(1 iPeeblesI 1200 ll : 
iHovle fc Vogelevl l200aT or the still controversial large- 
scale features feroadhurst et aljll99Ct iPercival et aljl200ll : 
iFrith et alJl200a7 justify the development of more complex, 
non-gaussian scenarios where linear, two-point predictions 
are in agreement with observations of the CMB, with the 
hope that simulations of galaxy formation in different envi- 
ronments and of the Lyman-a forest will agree with observa- 
tions at least a t the l e vel of ACDM dKauffmann et alJll999l: 
[ Springel etafl 1200 J: iMathis fc White! l2002t ICroft etall 
bool see also Ivan den Bosch et alJl2005t . 

We consider one such alternative to the concordance 
model, where primordial bubbles of true vacuum that 
formed in a first-order phase transition during inflation can 
survive to the present day and result in cosmological voids 
lO 119911: kiddle fc Wands! Il99ll : lOcchionero et al] Il997t 
loTchion^ro^^Arn^n^okl^94h ~This model is typically im- 
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plemented in the cont ext of so-called 'extended inflation' 
<La fc Steinhardtll989l). With reasonable values for the dis- 
trib ution of void s, iGjjflffhsetjflJ 12QQ-3, hereafter G03) (see 
also iBacciealupi et alJll997l : iBacciealupi fc Perrottal l200(jl 
have shown that the angular power spectrum of CMB tem- 
perature anisotropics in this primordial void model with 
ACDM-type cosmological parameters fits the observations 
at I <. 800. The signature of voids emerges at I J> 1000 
and, if as = 0.9, could account for the excess power seen 
by the Cosmic Background Imager (CBI) on cluster scales 
(I ~ 2500) if it is not due to an underestimate of the 
Sunyaev-Zel'dovich (SZ) (see fig. 3 of G03). While viable in 
ter ms of the tempera t ure a ngular power spectrum, we note 
that lCorosaniti et alJ (1200 J) derive constraints on the under- 
density and volume fraction of voids at recombination from 
the C OBE-D MR three-point correlation function and that 
IBacciealupi fc Perrottal <l2000l) predict non-gaussian features 
in the CMB temperature fluctuations that have yet to be 
tested against high-resolution maps. On galaxy scales, G03 
note that their void model may naturally account for the 
large underdensities seen in the local galaxy distribution 
which may be diffic ult to explain within the gau ssian ACDM 
paradigm ljPeeblesj200l|; lGottloeber et al. 2003, however see 
iMathis fc White! I2002T) . and use the observed typical local 
void radii and filling fractions to set two of the free param- 
eters of their void model. 

To assess this void model from another direction which 
could lead to simpler observational constraints, we simulate 
large-scale structure formation using collisionless dark mat- 
ter simulations of a network of compensated voids embedded 
in a ACDM cosmology. Due to the strong non-gaussianity, 
using concordance values for SIq , as results in a present- 
day cluster mass function (MF) that departs from the gaus- 
sian case. Non-linear structures develop much earlier in the 
void model, the result both of gravitational instability in 
the compensating dark matter shells surrounding the voids 
and of the large-scale motions triggered by the nonlinear 
evolution of the voids + shells systems. Integrating the clus- 
ter evolution up to z ~ 5, we show that simple number 
counts of SZ sources and optical depth to strong lensing 
are enhanced with respect to the ACDM scenario and that 
such observations may easi l y rule out the void model. Note 
that lAmendola fc Borganl lll994l) have compared the angu- 
lar two-point correlation function and the scaling of higher- 
order moments of a series of models with primordial voids to 
available observations of galaxy clustering. Their approach 
was however semi-analytical and lacked dynamical evolution 
of the structures. 

This paper is organised as follows. In Section 2, we recall 
the parameters of the void distribution and the approxima- 
tions we make, the same as G03. In Section 3, we discuss 
the set-up of the initial conditions, make simple checks and 
present results at z = 0. We deal with the high-z cluster evo- 
lution in Section 4, and motivate our choice of observables. 
We conclude in Section 5. 



2 PHENOMENOLOGICAL MODEL 
2.1 Parameters 

We take the parameters of the fiducial void model of G03. 
The background is a flat, dark energy dominated cosmology 



with h=0.7, Qo = 0.3, Ao = 0.7 and we normalize the am- 
plitude of the mass fluctuations so that as = 0.9. This set 
of parameters is close to the best-fit obtained from the com- 
bined WMAP + SDSS data llTeemark et al.ll2003rj| . We as- 
sume that fluctuations of the field driving inflation produce 
the usual gaussian adiabatic scale invariant perturbations, 
filtered as they reenter the horizon by a CDM transfer func- 
tion. Variants of extended inflation predict that the first 
bubbles to nucleate during the phase transition can reach 
cosmological sizes at recombination. The cumulative num- 
ber density of the voids is taken to be : 

N v (>r) =Ar~ a (1) 

where r is the physical void radius, A a normalisation con- 
stant adjusted to match the present-day filling fraction of 
voids seen in galaxy redshift surveys, and the exponent a 
is related to the gravitational coupling u) of the inflaton if 
extended inflation can be described with a Brans-Dicke for- 
mulation : 



Solar system experiments require uj > 3500 (Will 2001) and 
we will take here a — 3. The lower and upper cutoffs in 
the void radii r m j n . max are c hosen to ag ree with red shift 
survey s. IPlionis fc Basilakosl J2002T) and iHovle fc Voeelevl 
( 2002) measured the typical size of voids in the PSCz and 
UZC surveys and found that up to half of the volume of 
the universe is underdense with <5p ga i/p g ai <J —0.9. The 
void radii they obtain range from 10 to 30 h~ l Mpc with 
an average of ~ 15 ft -1 Mpc. Following G03 we assume 
r-min, max = 10, 25 h' 1 Mpc. (The contribution to the CMB 
temperature anisotropy angular power spectrum on CBI 
scales is mostly due to the voids with the largest radii; an 
analysis constraining r max using the CMB data alone is in 
progress.) Examples of physical motivations for the lower 
and upper cutoffs are that (1) on subhorizon scales, matter 
flows relativistically back into the voids after inflation dur- 
ing radiation domination, and suppresses the growth of small 
voids and that (2) the tunnelling probability of inflationary 
bubbles is modulated through the co upling to another field, 
resulting in a maximum radius. (See lOcchionero et al.lll99^ 
for a model of extended inflation involving two scalar fields 
and which gives a lower cutoff.) Last, we take an observed 
z = void filling fraction of /voids = 40%. 

2.2 Approximations 

In this section, we briefly review the assumed profiles and 
spatial distribution of the voids. 

We neglect the contribution of baryons as our purpose 
is to obtain statistics for massive clusters, where the dynam- 
ics is controlled by the dark matter. As the CDM becomes 
non relativistic early in the expansion, it is expected to travel 
only minimally into the voids which reenter the horizon later 
and are of interest for structure formation. On the other 
hand, due to the tight coupling to photons, baryons will fill 
in voids which are within t he horizon before decou pling at 
the adiabatic sound speed tLiddle fc Wands! Il99ll) . before 
gravitation takes over from radiation pressure. If at some 
point the mass of baryons inside the voids reaches a sub- 
stantial fraction of the mass of the shell, the growth of the 
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void radius will be slower than expected for a fully empty, 
compensated region. We simply mention that the precise dy- 
namics of void filling by the baryons is complex and needs 
further study when comparing for instance the high-redshift 
distribution of voids with the clustering of the Lyman-a for- 
est, and that the precise density profile of the voids depends 
on the physics of reheating and of the subsequent filling by 
baryons. 

Smooth void density profiles have been considered in an- 
alytical work: rounded step functions (iHoifman et alJll983t 
i Martel fc Wassermanl |l990) , 

exponentials lIHausman et alJ 
1983) or periodic fun ctions ( Baccigalupi 1998). However, 
with the exception of iReeos fc Gelled <ll991f) who also use 
a smooth initial profile, numerical simulations of void evolu- 
tion using pa rticles often approximate voids as "top-hat" un- 
derdensities JPubinski et al.lll993l: IWhite fc Ostrikerlll9961 
iRobinson fc Baked l200d hereafter W90, R00). These au- 
thors do not consider any particular physical model, but 
rather use generic templates of single voids or void networks. 
Here, to follow G03 we suppose that the voids at decoupling 
are spherical "top-hat" underdensities with <5pdm/pdm = 
— 1 surrounded by a thin compensating shell of dark matter. 

In an EdS cosmology, the shell density profile can be 
exactly derived from the self-simi lar solution of the evolu- 
tion of a spherical underdensity flBertschingerl ll98iL here- 
after B85). We assume that all the matter swept up during 
the expansion of the void ends up in the compensating shell 
around it. This behaviour results naturally from the expan- 
sion of an underdense region in an EdS cosmology (where 
the thickness of the shell is very slowly growing, see B85), 
a very good approximation to ACDM when we start the 
simulations. Recall that G03 assume for simplicity an EdS 
cosmology to compute the void contribution to the angular 
power spectrum of the CMB anisotropics, and add the re- 
sult to the power spectrum of the concordance ACDM model 
employed here. This makes sense, because (1) most of this 
contribution comes from voids close to the last scattering 
surface (LSS) where ACDM is similar to EdS and (2) the 
contribution to the power spectrum varies smoothly with 
I, reducing the small impact of the correction for angular 
diameter distance even further. 

As in G03 we ensure that the voids do not initially over- 
lap, although there seems to be no physical motivation for 
such a restriction. In practice, given the starting redshift 
we choose, our results depend very weakly on this hypothe- 
sis. Finally, the positions of centers of the voids are initially 
uncorrelated. 



3 SIMULATING THE VOID NETWORK 
3.1 Initial conditions 

We focus on two collisionless simulations of side 200 
h' 1 Mpc using 128 3 particles, carried out with the pub- 
licly available N-body tree-SPH code GADGET without hy- 
drodynamics . The first simulation, called G, is a gaussian 
ACDM model with the above parameters, the second, V, 
is the ACDM + voids non-gaussian fiducial model of G03. 
Except for the total initial displacement field, all simulation 
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parameters are similar in the two cases. The simulations em- 
ploy a Plummer softening length e = 0.08 h~ x Mpc which 
was kept fixed throughout in comoving coordinates. 

The size of the box is a compromise between the neces- 
sity of having enough primordial voids covering the whole 
radius range at z = 0, [10 25] h" 1 Mpc, and the mass res- 
olution. (The comoving radius range is [3 6.3] ft -1 Mpc at 
the starting redshift Zi n j t .) Even independently of the halo 
resolving power, A r part s = 128 3 is a stringent lower limit to 
the number of particles as simulations need to (1) propagate 
information about the smallest voids present at the starting 
redshift, (2) have sufficient initial power in the cosmologi- 
cal displacement field with respect to the power due to shot 
noise as one approaches the Nyquist frequency of the parti- 
cles not to alter significantly the formation of the smallest 
haloes one can resolve and (3) ensure that 2-body scattering 
effects are not important. While condition (1) is easy to ver- 
ify with our number of particles, (2) and (3) are a little more 
complex, as these issues can be amplified because we use 
non-gaussian initial conditions for V. In addition, because 
of the strong primordial non-gaussianity, the formation of 
structure may also depend on the type of the initial dis- 
tribution of particles which is employed before applying the 
initial displacement field, i.e. whet her one starts from a grid, 
a ran dom (Poissonian) , or a glass ( Bau gh et al ]|l995UWhitel 
1996) distribution. In the case of a grid, shot noise is mini- 
mal but the mesh introduces a characteristic scale and three 
preferred directions. A fully random, Poissonian distribution 
has no preferred direction but shot noise is significant. As 
a compromise between the two, we have used a glass dis- 
tribution, which does not have any preferred direction, but 
with intrinsic power spectrum rising as P(k) oc fc 4 . We find 
haloes with a friends-of-friends algorithm jDavis et alll985l) 
with a z — linking length bl — 0.164 times the mean in- 
terparticle separation I. At redshifts z — 1, 2 and 3, we have 
used the usual EdS linking length parameter 6 = 0.2 since 
ACDM behaves as an EdS cosmology at these epochs (this 
corresponds we ll to the 6 oc (Ac/O)" 1 ^ 3 scali ng proposed by 
lEke et al.lfl99rl see also I.Tenkins et aLlEoOlh . We keep only 
the groups with more than 10 particles and the minimum to- 
tal halo mass we can resolve is M min = 3.16 x 10 12 /i _1 M Q . 
To address the issue of the choice of the initial distribu- 
tion, we have checked that using a grid rather than glass for 
the initial distribution of particles in V does not change the 
z — MF of dark matter haloes. To deal with the impact 
of the level of shot noise compared to the amplitude of the 
initial cosmological perturbations, we verify below that the 
simulated z — MF of G matches analytical results (this 
is sufficient as we will show that the initial power spectrum 
of G is smaller than that of V over all simulated scales). 
Finally, we have found the z = number of dark matter 
haloes of V to depend on the starting redshift when using 
only Arparts = 64 3 particles, a signature of unphysical res- 
olution effects, but to have converged with Aparts = 128 3 
particles. 

In G, the initial displacement field d is given by the 
usual Zel'dovich approximation applied to the CDM power 
spectrum, normalised to a present as — 0.9. We compute d 
on a 128 mesh. In V, it is the same d for particles outside 
any primordial void and the displacement predicted by the 
similarity solution of B85 taken in the EdS regime (r oc t 4//o 
where r is the physical void radius and t the time) for all 
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particles that fall within a primordial void. Given the as- 
sumptions of paragraph 12.21 we put all particles falling in a 
void at its radius, and assign them the radial velocity of the 
expanding shell, following W90. In that sense, our V simu- 
lation would be close to that of a "spontaneous creation" of 
voids at Zmit- 

To follow G03, we use the similarity solution of an EdS 
universe (in particular the EdS time elapsed from Zimt to 
z = 0) to compute the initial radii and shell velocities of the 
voids, although we run the simulations with a ACDM back- 
ground. In doing so, we neglect differences expected with 
respect to EdS in the initial radii if voids were scaled back 
from z = using a theoretical solution for a ACDM cosmol- 
ogy. However, we verify below that this enables simulations 
to correctly reproduce single 25 /i" 1 Mpc radius voids at 
z — 0. Furthermore, the similarity scaling between void ra- 
dius and shell velocity remains strictly valid at zjnjt if it is 
high enough as EdS is then a very good approximation to 
ACDM (see B85). 

The starting redshift for simulations of cosmological 
models with gaussian initial conditions is usually set by re- 
quiring the maximum particle displacement to be ~ 30 per- 
cent of the mean interparticle distance. An upper limit on 
Zinit results from condition (1) on ./Vparts and translates the 
lower limit on the amplitude of the power spectrum of the 
initial perturbations. To choose Zinit for their simulations 
of a non-gaussian model with primordial voids, R00 ensure 
that A 2 < 0.15 at the Nyquist frequency of the particles 
and that no shell crossing occurs when initially displacing 
the particles. This is feasible as they compute the displace- 
ment for all their particles only from the Zel'dovich scheme 
applied to a linear density field. Their density field is the su- 
perposition of a gaussian field and a distribution of mildly 
underdense spherical regions with <5 V oid ~ 0.1. Because we 
start our simulations with compensated empty voids, we im- 
mediately probe the non-linear regime on void scales. As a 
consequence, the MF of haloes that we obtain at z = 
in the void model could depend on Zi n i t : it is possible that 
some dark matter haloes that our simulation can resolve 
have formed before Zj n j t by fragmentation of the void shells, 
and starting at zjnit could ignore them. Therefore we con- 
servatively start the simulations soon after decoupling with 
a high Zi n it = 1000, still satisfying condition (1). We have 
checked that the z — halo mass function in V obtained 
from starting at Zimt = 1000 is similar to that obtained 
from starting at Zi n i t = 3000 (this is not true for a simula- 
tion of V with only A r parts = 64 3 particles). Independently 
of the shape of the one-point PDF of the initial overden- 
sity field, a worry from such a high starting redshift is that 
numerical integration of the equations of motion might sup- 
pr ess the growth of s mall-scale modes (see the introduction 
of IScoccimarrolll998f) . Again, we use below the z = halo 
MF of G to check that this does not affect our conclusions. 

The initial gravitational potential <E>tot ( x*) that we con- 
struct in the simulations can be linked to CMB observa- 
tions. On void scales in the simulation V, like the overden- 
sity field, $ t ot will deviate f rom a realisation of a gaussian 
random field. As an example, Ko matsu et al.l a2003i) look for 
a signature of non-gaussianity in the WMAP CMB tempera- 
ture fluctuations maps using variants of the bispectrum and 
Minkowski functionals. They then constrain the relative con- 
tribution of second-order deviations from a gaussian field 



expected in single-field inflation models (e.g. iGangui et alJ 
11994 IVerde et al]l2000ft to the standard deviation of the 
gravitational potential but on a scale larger than the typical 
size the voids studied here would have on the LSS. Never- 
theless, analysis of non-gaussianity on smaller scales will be 
possible soon with interferometric maps, and further con- 
straints on the fiducial ACDM + voids model will result 
from comparisons with the theoretical bispectrum obtained 
in simulated maps of CMB temperature fluctuations. 

Here we adopt a simpler approach to the non-gaussian 
contribution to the gravitational potential and estimate the 
relative amplitude of the gravitational potential energy of 
the gaussian component with respect to the total energy in 
the compensated voids. All energies are computed in proper 
coordinates; note that the total (kinetic and gravitational) 
energy of a compensated void, including the shell, is always 
positive (see B85). Recall also that in an EdS universe W is 
constant if the growth of the fluctuations is linear. We note 
E v the total energy of the network of voids and shells in 
the V simulation and W s the total gravitational potential 
energy of the gaussian density fluctuations as realised in the 
G (or V) simulation. At z; n it = 1000, we obtain analytically 
E v = 7.09 x 10 22 while W g = -7.79 x 10 22 ft" 1 M© km 2 s" 2 . 
In the remainder of this section, we will implicitly assume 
this unit when quoting energies. 

Our analytical estimates show that at the starting red- 
shift, the energy due to the voids is similar to that of grav- 
itational potential of the gaussian density perturbations. In 
addition, the probability distribution function of the gravi- 
tational potential of a network of compensated voids over a 
homogeneous background is not gaussian (it is strongly posi- 
tively skewed). As a re sult, we expect / 2 j in the expansion of 
iKomatsu et alJ <|2003T) to reach values of order unity on the 
angular scale of the voids, but a precise calculation also sep- 
arating kinetic and gravitational energy of the voids+shells 
network is left for future work. 

In practice, because of shot noise and of transients from 
the Zel'dovich approximation, it is not easy to obtain a pre- 
cise value for the potential energy at high redshift directly 
from the particle distribution in the simulations (irrespective 
of a possible primordial non-gaussianity). We have checked 
that at z ~ 10 the predicted and measured values of W s 
agree within 20 percent. 

At z = 0, using the EdS similarity solution for the 
evolution of the network of voids and shells in a ACDM 
background, we obtain E v = 4.45 x 10 22 , and assuming the 
ACDM linear growth of the gaussian density fluctuations to 
be valid up to the Nyquist frequency of the particles we get 
W g = —6.07 x 10 22 . With this hypothesis of linearity in our 
whole box, the energy of the voids and shells drops at late 
times to 70 percent of the potential energy. In other words, 
from Zi n it to z — the analytical E v decreases by more than 
the factor Ain,ACDM/Ain,Eds(2imt — > 0) < 1 describing the 
evolution of W s . 

Directly in the simulations, we measure the total gravi- 
tational potential W sirn = -2.22x 10 23 in V and -1.85x 10 23 
in G. The excess factor of 2-3 between measurements and 
analytical predictions for W s im in G is due to the non-linear 
evolution of the density field. The lower VK s im measured in 
V compared to G can be a consequence of the slightly larger 
number of haloes that we find in V at z — 0. On the basis 
of the predicted E v at z = 0, one would expect W s im to be 
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Figure 1. Final particle distribution in a ACDM background 
simulation (without the cosmological gaussian perturbations) of 
a 25 h~ 1 Mpc - radius void centered on a 100 ft" 1 Mpc box. The 
slice shown has thickness 10 ft" 1 Mpc. The radii of the inner and 
outer solid circles are respectively the comoving input at Z[ n n (6.3 
h- 1 Mpc) and 25 ft" 1 Mpc. 

higher in V than in G. This is not measured, however, be- 
cause as seen in G non-linear effects determine the late-time 
value of W s im- 

For completeness, the total kinetic energy measured in 
the simulations right after setting the initial conditions is 
K tot = 1.44 x 10 16 and 1.60 x 10 14 in the V and G simu- 
lations respectively (Ktot = 1.19 x 10 23 and 1.02 x 10 23 at 
z—0) . The peculiar velocities assigned to dark matter parti- 
cles constituting the void-compensating shells is responsible 
for the large excess of initial kinetic energy in V. At late 
times, this excess is washed out in the dominant contribu- 
tion of the particles of virialized haloes. 

3.2 Consistency checks 

To check the effect of periodic boundary conditions, we sim- 
ulate with GADGET the growth from zjnjt of a single void 
with final radius 25 ft" 1 Mpc centered on a 100 ft" 1 Mpc 
box without the gaussian part of the displacement field, in 
a EdS universe. We use the EdS similarity solution to scale 
back to Zinit- At z = 0, we find that the total mass enclosed 
inside 99 % (98 %) of the inner shell radius of 25 ft" 1 Mpc 
is 9 % (3%) of that expected in a similar volume with mean 
density. The same test is repeated with the scaling method 
described above for simulations in a ACDM background and 
we find a similar "leaking" mass fraction, see Fig. In both 
cases a large fraction this remaining mass is contributed by 
particles of clusters formed on the shell. Down to z — 0, the 
particle distribution remains smooth beyond a 27 ft" 1 Mpc 
radius. 

We verified that the z — halo MF of the G simu- 
la tion gives results in agreement with the fitting formula 
of Ijenkins et alJ i200ll hereafter J01) and that it is robust 
against changes in z- m n. (Using iVparts = 32 3 rather than 
-Mparts = 64 3 or 128 3 for the G simulation results in an excess 
of low mass haloes if Zinit = 1000, and using N part s = 64 3 
rather than 128 3 for the V simulation leads to a redshift 
dependence of the MF: a factor 3 more massive haloes than 
for iV P a rt s = 64 3 and z init = 100). 

The validity of using the Zel'dovich approximation 
when setting up the initial conditions is checked on large 
scales by the growth of the largest modes of the simulation, 
in agreement with linear theory down to z = 0. We have 
finally performed two other simulations of the void model, 
changing both the gaussian random perturbation field and 
the positions of the voids, and found results very similar to 
those for V. 

3.3 Results at z = 

Fig.|5|shows the projected (2D) density of two slices of side 
200 ft" 1 Mpc and thickness 20 ft" 1 Mpc cut at the same 
position through the V and G simulations (left and right 
panels respectively) . Recall that the two simulations use the 



same initial gaussian displacement field, which in the case of 

V is combined with the displacement due to the primordial 
voids. The grey scale is the same in the two cases. Note the 
voids in the slice through V, apparent immediately below the 
centre of the picture and at the middle of the upper frame, 
together with a void network developing at the centre left. 

Fig. |3] gives the initial and z = real-space overden- 
sity power spectra (dotted and solid lines for the G and V 
simulations respectively). Also shown are the power spectra 
of the unperturbed glass initial particle distribution (dash- 
triple dotted line), and of the initial conditions for the G 
and V simulations (dashed and dash-dotted lines). All spec- 
tra have been divided by the linear growth factor for clar- 
ity. The diamonds with associated error bar s show the real- 
space g alaxy power spectrum as measured bv lTegmark et all 
l|2003al) from the DR1 of the SDSS survey. 

The imprint of the voids is significant at Zinit at k > 
0.1ft Mpc" 1 , as is clear when comparing the dash-dotted 
and dashed lines. At z = 0, non-linear power is larger in 
the void model than in the ACDM model, and departure 
from the evolved ACDM power spectrum occurs as early as 
k = 0.1 ft Mpc" 1 . Of course a precise comparison between 
the V power spectrum and the data would need to correct 
for bias, a quantity that might show a different behaviour in 
the ACDM + voids model than in ACDM. If a precise as- 
sessment of galaxy bias falls beyond the scope of this work, 
we note that at k = 0.25 ft Mpc" 1 , the dark matter bias 
of the ACDM + voids model with respect to the ACDM 
model is not more than 1.15. At the same scale, the bias of 
the evolved dark matter density field of the ACDM + voids 
model with respect to the galaxy distribution reaches 1.3, a 
high but plausible value. 

The main panel of Fig. 2] compares the halo MF mea- 
sured at z = 3, 2, 1 and from left to right in the G and 

V simulations (dotted and solid lines respectively) to the 
data and to the fitting formula of J01 calculated for the pa- 
rameters of the gaussian ACDM model we have simulated 
(dash-dotted lines). While the G simulation agrees well with 
J01, there are systematically more haloes in the V simula- 
tion at z = 0, for mass thresholds M min < 5 x lO 14 /!" 1 M . 
The excess is a factor 1.5 in the abundance. For higher mass 
thresholds, the abundance is similar in the two models. This 
is more clearly seen in the insert in Fig. 0] which gives the 
differential halo MF measured around 2 x 10 14 ft _1 M®. The 
excess of the non-gaussian model is apparent at all mass 
scales smaller than 3 x 10 14 ft _1 Mq. For higher masses the 
two curves are identical. We find one halo with M to t > 
10 15 ft _1 Mq in both the G and V simulations. (Recall that 
the gaussian part of the initial overdensity field is normalised 
to t he same present - day as = 0.9.) 

IBahcall fc Cenl Jl993l) give N >M = 2 ± 1 x 10" 6 
ft 3 Mpc" 3 at M = 4 x 10 14 ft" 1 Mq. Even assuming their 
uncertainty factor of 1.3 in the mass of rich clusters, this 
normalisation falls significantly below both our V and G 
mass functions. The present-day MF of the V simulation 
is therefore only marginally consistent with this normalisa- 
tion, but the level of disagreement is the same as that of 
a ACDM model with concordance cosmological parameters 
and as = 0.9. Diamo nds with error bars o n Fig. Ulsh ow the 
more recent data of IBahcall et al.| (120031 ): fBahcall fc Bode! 
(2003) who constrain the amplitude of mass fluctuations 
using the "optical" mass function of clusters selected from 
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Figure 2. Left and right panels: projected 3D density at z = in slices cut at the same position through the ACDM + 
voids and gaussian (ACDM) simulations respectively. The side and thickness are 200 and 20 h~ 1 Mpc. The gaussian initial 
displacement field is the same in both simulations. The grey scale is the same for the two panels. Note on the left panel the 
voids right below the centre of the picture, at the middle of the upper frame, and a void network developing a "honeycomb" 
structure at the centre left. 




Figure 3. Real-space matter overdensity power spectrum at z = of the gaussian (ACDM) and non- 
gaussian (ACDM + voids) simulations : dotted and solid lines respectively. The dashed and dash-dotted 
lines show the same respective quantities at z = 1000, after division by the linear growth factor for clarity. 
The dash-triple dotted line is the power spectrum of the raw glass file, for comparison. Note the signature of 
the voids in the initial conditions at k > 0.1 h Mpc -1 and the stronger z = non-l inear power in t he void 
model. Diamonds with error bars show the real-space galaxy power spectrum from lTegmark et alJ J2003al) 
(a more direct comparison between the ACDM + voids simulation and the data would need a detailed model 
for galaxy bias). 



th e SPSS EDR. They fi nd a normalisation similar to that 
of iBahcall fc Cenl (1993), with best- fit values as = 0.9 and 
Qo = 0.19. From a combinat i on of X-ray and optical observa - 
tions. I Vikhlinin et al.l j2003t) : IVoevodkin fc Vikhlinin! (120031) 
derive constraints on the cosmological parameters f2o, Ao 
from the evolution of the cluster baryon mass fraction up 
to z ~ 0.5, requiring as to be set by the observed am- 
plitude of the z = baryon mass function. Their z = 
and 0.4 < z < 0.8 total mass functions for massive clus- 
ters are reported with triangles and squares respectively on 
Fig. [I] (we have omitted their error bars for clarity, and cor- 
rected for a weak mass depende nce of the cluster baryon 
fr action. seelVikh linin e t al.ll2003T) . T he z = mass functio n 
of IVikhlinin et alJ J2003 ) agrees with lBahcall et al.1 (2003). 
At high redshift, the observed mass function has a lower 
normalisation than both the G and V z = 1 mass functions, 
a repeat of the z — differences. 

In fact, it is clear from Fig. [I] that the difference of the 
fiducial model of G03 with respect to ACDM is much more 
significant in the intermediate and high-z mass functions. 
Unfortunately, measurements of the cluster mass function 
at high redshift z > 1 from X-ray observations are sub- 
ject to uncertainties in, for instance, the redshift evolution 
of the L-T and M-T cluster scaling relations (see however 



IVoevodkin fc Vikhlininl 120031) . In the following section we 
will consequently use integrated SZ cluster counts and the 
abundance of giant arcs to probe the signature of the void 
model, rather than the cluster mass function. 

We conclude this section with two additional tests at 
z = 0, first changing the parameters of the void model and 
then comparing with two gaussian models. 

We have simulated another primordial void model V' of 
the same type as V but with a much steeper spectrum for 
the distribution of initial void radii (using i? m m,max = 10, 
40 h' 1 Mpc, a = 6, /voids = 40%). G03 show that this 
model also agrees reasonably well with high and low-Z mea- 
surements of the angular power spectrum of the CMB tem- 
perature fluctuations. The V' z = cluster MF shows a 
large excess compared to the V MF (an order of magnitude 
more massive clusters and as ~ 2) , so that this set of void 
parameters is directly ruled out. 

Second, it is necessary to check the impact of non- 
gaussianity by comparing to a model with gaussian initial 
conditions but with an initial power spectrum similar to that 
measured in the void model right after setting the initial 
conditions. For this purpose, we have simulated two gaus- 
sian models V\ and V% with the same cosmology and starting 
redshift as V and G but taking as initial power spectrum the 
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Figure 4. Main panel: mass functions of the gaussian (ACDM) and non-gaussian (ACDM + voids) simula- 
tions: dashed and solid lines respectively, measured at z = 3, 2, 1 and from left to right. The dash-dotted 
line is the fitting formula of J01 to ACDM. The diamonds with error bars show the "optical" mass function 
that lBahcall et all l2003h obtain using the EDR of the SDSS. The triangles and squares give the cluster total 
mass functions at z = and 0.4 < z < 0.8 respect i vely. They have been derived from the cluster baryon 
mass function constructed bv lVikhlinin et a! ] i2003ft : IVoevodkin fc Vikhlininl <2003h using a combination of 
optical and X-ray data (their error bars are omitted for clarity) . Insert : differential mass functions dN /d In M 
of the gaussian (ACDM) and non-gaussian (ACDM + voids) simulations: dashed and solid lines respectively, 
measured at z = up to 5 X 10 14 h _1 Mq. 



dash-dotted line -P V oids,mit of Fig. [3] Vi and Vi correspond 
to two different normalisations of the input power spectrum: 
for Vi we have taken P V oids,lnit of Fig. [3] divided by the lin- 
ear growth from z = 1000 to z — 0, while for Vi we have 
normalised so that the simulation has as = 0.9 at z = (the 
adopted normalisation for Vi results in a% ~ 4 at z — 0). We 
have found the z — cluster mass function of Vi to largely 
exceed that of V and G (there is a factor of 5 more haloes 
with masses M to t > 4 x 10 14 fr _1 Mq) and that of Vi to be 
abruptly cut at masses Mtot ~ 10 14 /i _1 Mq. In addition, 
the matter power spectra measured at z = in Vi and Vi 
retained the strong feature ( "bump" ) seen in the initial con- 
ditions at fcvoid > O.lh Mpc^ 1 as a step-like increase at fc vo id 
which was only little modified by the late-time non-linear 
evolution. This is clearly ruled out by the data. To sum- 
marise, the non-gaussianity of the primordial void model is 
necessary to approximately reproduce the z = observed 
cluster mass functions and mass power spectra: employing 
the initial power spectrum of the void model in a gaussian 
primordial density field results in a large mismatch to the 
data. The large peculiar velocities of the compensating shells 
surrounding the voids and associated to the scale fcvoid are 
not realised in a gaussian model with same initial power 
spectrum as V. As a result, the pattern at fcvoid in the power 
spectrum is more stable and persists longer in our Vi and V2 
tests than in V , where large velocities may dilute/broaden 
the feature. The particular shape of the power spectrum 



resulting from the void network is not transposable to a 
gaussian initial probability function. 



4 DERIVING OBSERVATIONAL 
CONSTRAINTS 

In this Section, we make predictions for the thermal SZ ef- 
fect and for simple statistics of strong lensing. We show that 
they differ substantially in the primordial void model from 
their values in a ACDM cosmology. We note here that other 
observables like the cosmic shear or the clustering of the 
Lyman-Q forest could also bring out the presence of primor- 
dial voids, but they are more complex than the former and 
may be affected by biases due to the non-linear evolution of 
the power spectrum. 

4.1 SZ source counts 

We first estimate the counts expected from th e detection 
of the cluster thermal SZ effect, up to z = 5. iKav et ail 
teOOll) make detailed analytical predictions for the SZ num- 
ber counts expected for the Planck satellite using large sim- 
ulations of cluster formation in gaussian CDM cosmologies. 
We follow the same approach but with simplifying assump- 
tions. We use 30 simulation outputs for G and V from z = 
to z ~ 5. Each dump has a comoving size 200 hT 1 Mpc and 
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Figure 5. 143 and 353 GHz resolution limit for the cluster ther- 
mal SZ effect in the simulations as a function of reds hift. The sen- 
sitivi ty of Planck to the SZ effect is of order 30 mjy IBartelmanrJ 
1200 ill , that of the planned Bonn-Berkeley APEX - SZ survey will 
reach 3 mjy at 150 GHz. 



we find their DM halos above the minimum mass threshold 
M m i n . For simplicity we assu me for each halo an i sothermal 
profile for the gas (see, e.g., iBarbosa et aljll996fl with to- 
tal mass Mgas = /b x Mtot where /b = 0.13 is the cosmic 
baryon fraction. Taking the clusters to be point sources, we 
compute the magnitude of the flux change observed against 
the CMB at the Planck satellite frequencies 143 and 353 
GHz, on both sides of the zero-point of the SZ thermal ef- 
fect (217 GHz). Fig. |5]shows the minimum SZ flux that the 
simulations can resolve as a function of z. In the following, 
we show results for fluxes S v > 10 mjy. 

Fig. |S| gives the expected redshift distribution of SZ 
sources with S v > 10 mjy in the ACDM + voids and ACDM 
models (solid and dotted, dashed and dash-dotted lines re- 
spectively), for 143 and 353 GHz. The curves have been 
fitted by polynomials in the relevant range. While there are 
virtually no sources with S143, 353 ghz > 10 mjy at z > 1.6 
(z J> 2.1) in G, the distribution of such objects in the void 
model is different. At 143 GHz, there are more sources in 

V at low redshift compared to G (a factor 1.3 in excess at 
z = 0, reaching more than 2 at z = 1), and the distribu- 
tion of sources extends to z — 2.6 on the window shown. At 
z > 4.3, beyond the "dip", the differential number counts 
of sources in V at 143 GHz again exceed 0.1, and keep in- 
creasing to z ~ 5 and earlier. At 353 GHz, SZ sources are 
intrinsically brighter and the differential number counts in 

V exceed 0.1 all the way to z ~ 5: the number counts de- 
crease from z = to z ~ 3.4, then increase to z = 5 and 
beyond. At z = 0, there is an excess by a factor 1.6 over 
the same V counts at 143 GHz, up to z ~ 2 where the ra- 
tio increases. At z = 5 the ratio between the 353 and 143 
GHz differential counts in V is 2.5. With respect to G, the 
353 GHz differential number counts in V are higher by 60 
percent at z — 0, and by more than an order of magnitude 
at z = 1.8. At low redshift z ~ 0, the excess of sources in 

V compared to G is accounted for by the excess of clusters 
of mass Mtot ~ 10 14 /i _1 M© seen Fig. The slower de- 
crease in the number counts between z = and 2 = 3 in V 
compared to G is the direct consequence of the slower late- 



Figure 6. Redshift distribution of the sources with flux 
S143 353 GHz > 10 mjy in the ACDM (dashed and dash-dotted 
lines) and ACDM + voids (solid and dotted lines) simulations. 



time evolution of the cluster mass function in V compared 
to G. Finally, the increase in the differential number counts 
at z > 4 in V is due to the flattening and decrease of the 
angular diameter distance at such redshifts in ACDM, to- 
gether with the presence of massive, sufficiently hot haloes 
at these epochs in the void model. 

In reality, there will be a residual number of sources 
with S143, 353 GHz > 10 mjy at z > 5 in the void model, 
as the results shown here provide only a lower limits to the 
counts, but the difference with respect to ACDM is already a 
factor of 2 at z — 5 and SZ surveys may therefore falsify the 
void model. The detection of sufficiently bright high-redshift 
SZ sources, on the other hand, would be a hint towards non- 
gaussianity. 

To see how this conclusion varies with sensitivity, Fig.|7| 
gives the number counts of SZ sources expected (up to z = 5) 
in the ACDM (dashed and dash-dotted lines) and ACDM + 
voids (solid and dotted lines) models as a function of the 
flux threshold. Above a flux limit of 10 mjy, there are ~ 2 
(~ 2.5) times more integrated counts at 143 (353) GHz in the 
ACDM + voids model compared to ACDM. T he sensitivity 
of Pl anck to the SZ effect is of order 30 mjy iBartelmanrJ 
2001) as shown by the vertical dash-triple dotted line, and 
the corresponding counts are enhanced by a factor 1.8 (2) 
at the two frequencies. For sources brighter than 100 mjy 
however, the excess in the number counts of sources in the 
primordial void model compared to ACDM drops to a fac- 
tor 1.5 (1.25) at 143 (353) GHz, a ratio maybe too small to 
falsify the void model against ACDM with upcoming obser- 
vations, even if a more detailed computation is needed at 
this level. Among many ground-based examples, the Bonn- 
Berkeley APEX - SZ survey 2 will cover 100 square degrees 
and is expected to reach a sensitivity of 3 mjy at 150 GHz. 
This 300-hour experiment will detect ~ 1000 clusters with 
mass Mtot > 2 x 10 14 /i _1 Mq up to z — 2 (assuming a 
ACDM model). It will provide strong constraints on non- 
gaussian models like primordial voids. 



'■ http : //bolo . berkeley . edu/apexsz/ 
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Figure 7. Number counts of SZ sources to z = 5 in the ACDM 
(dashed and dash-dotted lines) and ACDM + voids (solid and 
dotted lines) models as a function of the minimum flux S^min 
at 143 and 353 GHz. Note the enhancement (factor ~ 2.5) of 
the counts above S„ jm i n = 10 mJy at 353 GHz in the model with 
primordial voids. The vertical line shows the sensitivity of Planck; 
at this level, number counts differ by a factor ~ 2. 



4.2 Strong lensing 

The statistics of strong amplification of the images of back- 
ground sources by massive concentrated structures is a pow- 
erful probe of the cosm ological parameters (for example, 
iBartelmann et ai1ll99ct) . Clearly, they will also be directly 
affected by the primordial non-gaussianity considered here. 

We assume that only the clusters that we can resolve 
in the simulations, (A/tot M m i B ) produce strong lenses. 
We model the h aloes with singular isothermal spheres (SIS) 
JPeacockl 11983) and measure their velocity d ispers ions di- 
rectly from the simulations. IPerrotta et al.l <l2002l) derive 
analytical predictions with NFW profiles for the differen- 
tial amplification probability P(A) and compare them to 
the predictions obtained with SIS profiles. For sources at 
z s = 4 and 7, at A <; 3 they find their ACDM model with 
NFW profiles to be more efficient than with SIS profiles. At 
A > 7, SIS profiles yield the higher P(A), but the probabil- 
ities for NFW and SIS profiles stay within a factor 3 of each 
other even at larger amplifications. We note here that we 
also neglect the impact of substructure and asphericity on 
the cross section for strong lensing, which may be significant 
but n eeds better understanding (see, e.g.. iMeneehetti et alJ 
12002^ 

have computed the number of 



IBartelmann et all ( 199S 
giant arcs expected on the whole sky for a series of gaussian 
CDM models and shown that only open CDM models could 
reproduce the total number of arcs seen in the EMSS sam- 
ple. In particular, their open CDM model produced an order 
of magnitude more giant arcs than their ACDM model, as- 
suming a redshift z s = 1 for their sources. In a recent work, 
IWambseanss et all i2003l) study the effect of the redshift dis- 
tribution of sources in the predicted number counts of giant 
arcs in a ACDM cosmology. They show that the lensing op- 
tical depth is a very steep function of the source redshift for 
0.5 < z s < 2.5 and that it increases further for z s > 2.5. 

In fact, the order of magnitude discrepancy between the 
EMSS counts and the ACDM prediction is solved if one takes 



a more extended distribution of sources to high redshifts. 
This is more realistic as a large fraction of the giant arcs 
Idadders et alJ i2003t) count in the Red-Sequence Cluster 
Survey (RCS) have high redshifts: 1.7 < z alc < 4.9. Putting 
a third of the sources at z s — 1.5 rather than all sources at 
z s — 1 increases the ACDM predicted counts by a factor 3, 
and putting 7 percent of the sources at z s > 3 increases the 
counts by a facto r of 7, bringing ACDM in agreement with 
the observations JWambseanss et alJl20o3l . 

Because the lensing optical depth is such a steep func- 
tion of source redshift, one expects the large differences be- 
tween the ACDM and the primordial voids models in the 
mass function of clusters at z ~ 1 — 3 (probing high source 
redshifts) to significantly affect the number counts of giant 
arcs, even without classifying the counts along arc (source) 
redshift. In the following, for simplicity, we compute our 
strong lensing statistics putting all our sources at high red- 
shifts: z s = 3 and then at z s = 5. Detailed comparison to 
observations would require a more realistic model for the dis- 
tribution of source redshifts, but our plots show the level of 
discrepancy expected between ACDM and ACDM + voids. 

Fig. |H] shows the cumulative probability P(A > A m in) 
that a line of sight has an amplification larger than A m i n in 
the ACDM and ACDM + voids cosmologies and for the two 
values of z s . We take the strong lensing regime to be A > 2. 
At Amin > 10, the cumulative probability for strong lensing 
is a factor ~ 2 (~ 4) higher in the void model than in the 
ACDM case, for z s — 3 (5). For sources at z s = 5, we note 
that the enhancement is also larger than the possible bias 
due to our choice of an SIS rather than NFW halo profile. 
We conclude that the optical depth to strong lensing is in- 
creased in the ACDM + voids model compared to ACDM, 
and that the number counts of giant arcs may differ by up 
to a factor of 4, hence overpredicting the observations com- 
pared to ACDM. The optical depth to high-redshift galaxies 
will be increased by a factor of 4 in the void model, and the 
the number counts of high-redshift giant arcs (with z arc ~ 5) 
will be increased by more than a factor of 4 in the ACDM 
+ voids model compared to the ACDM model, because the 
comoving density of bright (M* > 10 10 h~ x Mq) galaxies 
is also higher at z ~ 5 in the ACDM + voids model than 
in the ACDM model. (We verified this last point using a 
semi-analyti cal model for galax y formation similar to that 
described in lMathis et al][2002l) . If the density of the lensed 
population evolves similarly in the ACDM and the ACDM 
+ voids models, then the mean redshift of the lenses will be 
shifted to higher values in V compared to G: Fig. [5] shows 
the redshift distribution of the total optica l depth dr/dz t o 
strong lensing as defined in section 4.2 of |Peaco"ckl il999t) . 
for z s = 3 and 5. The total optical depth r is larger in the 
V than in the G simulation (note the factors ~ 2 and ~ 2.5 
decrease between the V and G curves at z s = 3 and 5 in 
Fig.®. 

Finally, we have assumed initially compensated voids 
surrounded by a thin shell growing as the underdensity 
expands in comoving coordinates. This large-scale config- 
uration could per se c onstitute an efficient lens. However, 
lAmendola et al.l il999t) show that only voids with radius 
larger than ~ 100 h" 1 Mpc today induce weak gravita- 
tional lensing with a signal to noise ratio greater than unity 
in observables like color-dependent density magnification or 
aperture densitometry. Even for strong underdensities like 
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Figure 8. Probability P(A) that a line of sight is magnified by 
A > A min for the ACDM (dashed and dash-dotted) and ACDM + 
voids model (solid and dotted) for two source redshifts (z B = 3 and 
5 respectively). We suppose that the only lenses are the massive 
DM haloes with M to t > A/ min = 3.16 X 10 12 /i _1 M , modelled 
as singular isothermal spheres. 
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Figure 9. Redshift distribution of the optical depth to strong 
lensing dr/dz contributed by our resolved haloes. The legend is 
the same as for Fig. Note the ratio of ~ 2 to 2.5 between the 
primordial void model and the gaussian ACDM cosmology. 



<5void ~ — 1, none of these weak lensing techniques will be 
able to find a significant signature of the voids considered 
here. 



5 CONCLUSIONS 

We have simulated cluster formation in a physically plau- 
sible non-gaussian primordial void model where empty and 
fully compensated bubbles surviving from inflation together 
with the gaussian adiabatic CDM-type perturbations pro- 
vide the seeds for the development of structure. This model 
shares the cosmological parameters of ACDM and possesses 
gaussian statistics on large scales while non-gaussianity only 
affects the one-point distribution on cluster scales. It is an 
attractive alternative to ACDM as it may explain the ex- 
cess of CMB temperature anisotropy power at I ~ 2500 



recently observed by CBI and as it could account for the 
large voids seen in the nearby galaxy surveys JCriffiths et alJ 
2003). While analysis of high-resolution CMB maps using 
higher-order moments will provide further constraints on 
primordial non-gaussianity, we have shown the evolution of 
the cluster mass function to be a strong constraint at low 
and intermediate redshifts. 

At startup, the total energy of the network of voids 
and shells in our simulation is comparable to the potential 
energy of the gaussian fluctuations. Even if strong under- 
densities are present on Mpc scales very early on compared 
to ACDM, we have found that our 200 ft -1 Mpc-side, 128 3 - 
particle simulations which start shortly after recombination 
can provide a reliable estimate for the cluster mass function 
of the non-gaussian model. 

The power spectrum of the void model measured after 
setting the initial conditions shows a strong feature char- 
acteristic of the voids at k J> 0.1 ft Mpc , which is then 
erased during the evolution by the non-gaussian initial con- 
ditions. The z — matter power spectra of the ACDM + 
void model is close to that of the gaussian model, with little 
additional small-scale power at k > 0.5 ft Mpc -1 . The z = 
mass function gives a similar number of massive clusters 
(Mtot > 4 x 10 14 ft _1 M Q ) in the ACDM + voids model 
and in the concordance ACDM scenario. Our high clus- 
ter abundance compared to observations is a consequence 
of the WMAP+SDSS normalisation. The evolution of the 
cluster mass function up to z ~ 3 in the void scenario dif- 
fers strongly from that of the ACDM model and is much 
more "efficient" than the sole z — mass function in distin- 
guishing between the two. The subs tantial evolution in the 
cluster baryon mass function seen bv lVikhlinin et all <l2003l) 
between z = and z ~ 0.5 needs confirmation before we 
can falsify the non-gaussian model which has little late-time 
evolution. 

In fact, better constraints on the void model are ob- 
tained from two other statistics at z > 1.5. We have first 
shown that the integrated number counts of SZ sources is 
higher by a factor 2.5 (resp. 2) in the ACDM + voids model 
compared to the ACDM model, for 353 GHz flux greater 
than 10 mjy (resp. 30 mjy, the resolution of Planck). We 
have then used the optical depth to strong gravitational lens- 
ing as another possible discriminant between the primordial 
void model and ACDM. We have shown that the number 
counts of high-redshift, z > 3 (resp. 5) arcs is expected to 
be more abundant by a factor of 2 (resp. 4) in the non- 
gaussian scenario, for the same underlying lensed popula- 
tion. Because the optical depth to strong lensing is a steep 
function of redshift, we expect the total number of giant arcs 
observed to be increased by a factor > 2 in the ACDM + 
voids model compared to ACDM. As a result, the ACDM + 
void scenario overpredicts the number counts of giant arcs 
seen in the EMSS compared to the concordance ACDM. 

An additional contribution to the CMB power spectrum 
at high £ will also be generated by the thermal SZ contri- 
bution from unresolved clusters. The joint limits on a non- 
gaussian contribution from clusters to the power spectrum, 
including the ROS AT/WMAP cross-correlation constraints 
JPieeo et al.ll2003l) and the CBI and ACBAR experiments, 
will be presented elsewhere. 

Not only do such simple tests give the opportunity to 
rule out a particular set of models, they also more generally 



Constraining C MB- consistent primordial voids with cluster evolution 11 



bring out simple constraints that have to be satisfied (e.g. 
with numerical simulations) as one proposes non-gaussian 
alternatives to the current paradigm. 
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